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NONSYMMETRIC  BALLASTIC  RANGE,  HEIGHT,  TIME -OF -FLIGHT 
AND  OPTIMAL  FLIGHT  PATH  ANGLE  COMPUTATIONS 
WITH  PROGRAMS  FOR  A  HEWLETT-PACKARD  65  CALCULATOR 


I .   The  Introduction  and  Apology. 

The  range  of  a  ballistic  missile  over  the  surface  of  the 
Earth  is  a  topic  which  is  discussed  in  almost  every  text  on 
astronautics  (Reference  1,  for  example).   Unfortunately,  none 
of  the  available  texts  contained  a  needed  procedure  for  the 
exact  solution  of  the  unsymmetric  ballistic  problem  (defined 
in  Section  II) ;  all  procedures  found  were  approximate  and 
were  based  on  a  symmetric  approximation  [See  References  2  and 
3/  for  example] .   The  author  of  this  report  faced  the  dilemma 
of  either  (1)  taking  the  time  to  solve  the  problem  afresh,  or 
(2)  taking  an  inordinate  longer  time  to  do  a  thorough  search 
to  find  a  solution  which  might  not  be  the  one  desired.   The 
iormer  course  of  action  was  taken  and  all  due  apologies  are 
hereby  extended  to  those  that  have  published  more  elegant 
solutions . 

All  of  the  computational  methods  are  outlined  in  Section 
III.   Instructions  and  a  sample  problem  for  the  Hewlett- 
Packard  65  Calculator  (HP-65)  are  given  in  Section  IV.   This 
HP-65  programs  are  in  Section  V,  and  the  development  of  the 
3   .   procedure  is  in  the  Appendix.   Familiarity  with 
Newtonian  two-body  theory  is  assumed. 


II.   The  Problem. 

Basically,  the  problem  is  to  compute  the  range,  S,  and 
time-of -flight,  tf,  of  a  ballistic  missile  given  the  booster 
cutoff  at  height  h„   after  a  vertical  ascent,  a  cutoff 

£j 

velocity  vT ,  and  a  subsequent  flight  path  angle   3   (Fig.  1). 

Li 

Given  that  a  solution  to  the  stated  problem  is  available,  the 
next  natural  question  to  ask  is,  "What  flight  path  angle   3 

will  obtain  the  maximum  range   S    ?"   The  solution  is  fairly 

3    max  J 

simple  in  the  symmetric  ballistic  problem,  i.e.  when   h  =  0, 
and  called  symmetric  because  the  launch  point   L  and  target 
point   T  are  symmetric  with  the  line  of  apses  of  the 
Keplerian  orbit  (i.e.   <j>,  =  cj>2   in  Fig.  1).   The  unsymmetric 
ballistic  problem  is  that  in  which  h_  *  0 .   Here  we  assume 
that  h„   is  given,  the  Earth  with  radius   r_,   is  spherical, 
non-rotating,  and  has  no  atmosphere.   The  effect  of  the  Earth's 
rotation  may  be  calculated  figuring  how  far  the  target  has 
moved  during  the  time  of  flight  of  the  missile.   Explicitly, 
we  assume  that  Keplerian  orbit  is  obtained  at  height   h_  with 

Ei 

terminal  booster  cutoff  velocity  v_   and  flight  path  angle 
3 (or   3T  )  ,  and  the  target   T   is  at  the  Earth's  surface. 

Li 

The  computation  of  range  and  time-of -flight  for  given  h_, 

c 

r  ,  v  ,  and   3   can  easily  be  obtained  from  the  equations  of 

elliptical  orbits.   The  determination  of   3     is  somewhat 
r  opt 

complicated,  however,  and  requires  an  iterative  procedure.   In 
all  of  the  developments  we  assume  that   L  and   T   are  sepa- 
rated by  the  line  of  apses  (which  implies  that   0  <  3  ^  y  )  . 


Line  of  apses 


Unpowered  flight 


FIGURE  1 


III.   The  Computational  Procedures. 


A.   Range 


1.   Let   rL  and  rT  denote  the  distance 

of  the  launch  point   L   and  target  point 
T   from  the  center  of  the  Earth.   Usually 
we  let   rL  =  rE+hE  and  rT  =  rE 
although  r   and   r   need  not  be  so 

Li  1 

defined  in  general. 


2.   From  rT   and  vT ,  compute  the  semimajor 

Li  Li 


axis  of  the  ellipse, 


a  = 


ft)- 


2 

where   y  =  k   and  k   is  the  (Gaussian) 

gravitational  constant  for  the  Earth.   We 

have  assumed  that  the  mass  of  the  missile 

is  negligible  with  respect  to  the  mass  of 
•  the  Earth. 

3.   Using  the  flight  path  angle   3,  compute 
the  elliptical  semiparameter 


2 

P  =      (rLvLcos6)    /]i    . 


4.   Compute  the  eccentricity 


e  = 


5.   Compute  the  true  anomaly  of  the  launch 
point 


f-  =   cos 


'  [%  - ')] 


where   0  <  f   <  tt  .   Note  that  the  line 

Li 

of  apses  can  be  said  to  have  true  anomaly 

f  =  IT  . 

6.   Compute  the  true  anomaly  of  the  target  point 


f  =  2tt  -  cos 


-1 


[*fc  - ')] 


7.   The  surface  range  of  the  missile  is  then 


S  ■  rE(fT  "  fL>  • 


B.   Time-of -Flight 

1.   Compute  the  eccentric  anomaly  of  the 
launch  point, 


ET  =   cos 

Li 


-1 


U  -  '4 


where   0  <  E   <  tt 


2.   Compute  the  eccentric  anomaly  of  the 
target  point, 


E      =    2tt    -   cos 


-1 


M'  -  *)] 


where      tt    <   E_    <    2tt 


3.   Using  Keplers '  equation,  compute  the 
time-of -flight 


3/2 


(ET  -  EL)-e(sinET  -  sinEL) 


C.   Height  of  Trajectory  &  Circular  Velocity 

1.   The  maximum  height  of  the  trajectory  above 
the  Earth's  surface  for  any  h„,  vT ,  and 

t,  Li 

3   is  given  by 


h    =  h  =  a(l  +  e)  -  r_ 
max  E 


2.   Circular  Velocity 

The  velocity  required  for  a  circular 

orbit  at   r   is 
j-i 


/v* 


vc 


v   will  be  used  in  Section  III  D 
c 


D.   Optimum  Flight  Path  Angle  &  Maximum  Range 

For  a  specified  h_   and  vT ,  a  value  of 

£i  Li 

3  can  be  found  such  that 


Smax  '  S(6opt'  "  T  S(S)  • 


In  this  unconstrained  optimization  problem 

as 


we  set 


0        0  =0     and  solve 

6  ~  3opt 


for   3   .  •   The  relationship  between   S   and 
$   is  given  by  the  computational  formulas  for 
the  range  in  part  A.   Although  an  explicit 
formula  for   3   is  unobtainable,  the  following 
exact  relationship  is  useful: 


(1)   cos^3 


opt 


(2-a  )  J  r  J 


2x    20 

)cos  3 


opt 


2y/l     1  \1  .  , 
— £-/ —  -  — I  sin 


opt 


V 


-fA   cos23  ¥    -    2M~  ~   ~\ 
\rTJ  opt   vL^rL   rT) 


where   a  = 


v 


In  the  symmetric  ballistic  problem  rT  =  r„,  ; 
the  equation  for   3   .   then  reduces  to 


(2) 


cos23 


opt    2_a2  ' 


which  is  the  result  given  in  References  1,  2,    and  3 


In  the  unsymmetric  case  Equation  (1)  is  used  to 

find   3   a.   as  follows: 
opt 


1. 
2. 


Solve  (2)  for   3   .  . 

opt 

Substitute   3   .   into  the  right  hand  side 


of  (1)  to  obtain  a  new  value  for   3 


opt 
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3.   Average  the  last  two  values  of   &ODt  to 

obtain  a  new  value  of   3   .   and  return 

opt 

this  value  of   BODt  to  step  2. 


4.  It  has  been  found  that  five  iterations  are 
sufficient  to  obtain  a  reliable  value  for 

Sopt  ' 

5.  Using   3ODt»  compute  the  maximum  range 
Smax  =  S^opt*   usin9  procedure  A. 


IV.   The  Illustrative  Example, 

The  HP-65  program  consists  of  four  magnetic  cards  which 
are  used  as  follows: 

Card  1:   This  card  must  be  entered  first.   It  is  used 
to  set  the  two  constants   r„   and   y   for  any 
choice  of  four  distance  units  and  any  choice  of 
four  time  units.   The  distance  units  are 
kilometers,  meters,  statute  miles,  and  nautical 
miles;  the  distance  unit  is  entered  by  pressing 
key  A,  B,  C,  or  D,  respectively.   The  display 
shows  the  equatorial  Earth  radius   r^   in  the 
respective  unit.   The  time  units  are  seconds, 
minutes,  hours,  and  days;  the  time  unit  is 
entered  by  pressing  key  A,  B,  C,  or  D,  respec- 
tively.  The  display  shows  the  constant   y   in 

3        2 
the  appropriate  (distance)  /(time)   unit.   Note  - 

after  entering  program  Card  1,  the  distance  unit 
must  be  entered  before  the  time  unit  is  entered 
since  the  keys  A  through  D  perform  a  dual  func- 
tion.  The  basic  definitions  used  are: 

1  statute  mile  =  1.609344  kilometers (km) 
1  nautical  mile  =  1.852  km 
r„  =  6378.160  km 
y  =  398603  km3/sec2 


Card  2:   This  card  computes   S,  h.  ^  ,  and   v   given 

max         c 

h_,  vT   and   3.   It  must  be  used  after  the 

E      Li 


execution  of  Card  1  and  it  may  be  used  before 
and/or  after  the  execution  of  Cards  3  and  4.   If 
it  is  used  after  Card  4  to  compute  smax'  then 
$   .   as  computed  by  Card  4  need  not  be  re-entered 
Otherwise,  in  any  order:   enter   h„   and  press  A, 

IS 

enter  v_  and  press  B,  enter   6   and  press  C.   To 

Li 

find  the  range  S,  press  E.   Then  to  find   h.    , 

3      c  max 

press  R/S;  and  then  to  find  v  ,  press  R/S  again. 
Note:   If   h_,  vT ,  and   8   have  been  entered 

Ei      J-i 

using  Card  3 ,  then  they  do  not  have  to  be 
reentered  when  using  Card  2.   Also,  any  new  value 
of  h_,  vT ,  or   3   can  be  entered  and  the  range 

Ei      Li 

computed  without  having  to  re-enter  unchanged 
values. 

Card  3:   This  card  computes  the  time  of  flight,  tf  . 
It  must  be  used  after  the  execution  of  Card  1 
and  it  may  be  used  before  and/or  after  the 
execution  of  Cards  2  and  4.   If  it  is  used 
after  Card  4,  then   8   .   as  computed  by 
Card  4  need  not  be  re-entered;  if  it  is  used 
after  Card  2,  then  only  new  quantities  need  be 
entered.   Otherwise,  in  any  order:   enter   h_ 

Ei 

and  press  A,  enter  v_  and  press  B,  enter   8 
and  press  C.   To  find  the  time  of  flight,  press 
E. 


10 


Card  4:   This  card  computes   $   . .   It  must  be  used 

opt 

after  h   and  v   have  been  entered  from 

either  Card  2  or  Card  3.   No  entries  are 

required.   To  compute   3   t/  press  any  key  A 

through  E;  five  iterations  will  be  performed 

and  the  resulting  value  of   3     will  be  dis- 

opt 

played.   To  perform  five  more  iterations,  press 
R/S.   To  monitor  the  intermediate  iterations 
press  keys  f  SF2.   Each  intermediate  result  will 
be  displayed.   To  continue,  press  R/S.   To  dis- 
able the  monitoring,  press  keys  f   SF2. 


Using  the  four  program  cards  as  described  above  find 
the  range  and  time-of-f light  for  a  ballistic  missile  that 
has  a  velocity  of  10000  knots  at  a  cut-off  altitude  of  50 
nautical  miles  for  a  flight  path  angle  of  30   and  40  . 
Also,  find  the  maximum  possible  range  and  time-of-f light 
for  the  same  missile. 

Proceed  as  follows: 


Instruction 
Enter  Card  1 
Distance  unit  is  n.mi. 
Time  unit  is  hours 
Enter  Card  2 
Enter  cutoff  altitude 
Enter  cutoff  velocity 
Enter  flight  path  angle 
Compute  range 


Data 

Key 

'  D 

C 

50 

A 

10000 

B 

30 

C 

E 

Display 


rE=3443.  93n.ini 
y  =  8.13247xl0i:L 


r  =3493.93n.mi 


S=1928 .08n.mi 
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Step 

8a 

8b 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

Instruction 


max 


Optional:   Display  h 

Optional:   Display  v 

Enter  flight  path  angle 

Compute  range 

Enter  Card  3 

Compute  time-of-f light 

Note:   h_  -  50,  vT  =  10000,  and 
3  =  40°   since  no  further 
entries  have  been  made. 

Enter  flight  path  angle 

Compute  time-of-flight 

Enter  Card  4 

Compute   $   .   (5  iterations) 

Make  5  more  iterations  to  be  sure 

Enter  Card  2 

Compute   S„a„ 
max 

Enter  Card  3 

Compute  time-of-flight  for 
maximum  range. 


Data 


40 


30 


Key 

R/S 

R/S 

C 

E 

E 


Display 


*i   =344.74n.mi 

max 

v  =15256. 47knot 
c 


S=1956.60n.mi. 


tf=0.3043hrs. 


c 

E 

0.2504hrs. 

E 

8opt=36?0895 

R/S 

Bopt=36°°895 

E 

S_  =1973.95n.B 

max 

W'-0-2846 
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V.      The  Program  Listings. 
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CARD  1 


Set  units  of  distance  and  time 


display  r   and 


10 


Key 

Code 

Entry 

Shown 

LBL 

23 

B 

12 

EEX 

43 

3 

03 

CHS 

42 

GTO 

22 

0 

00 

LBL 

23 

C 

13 

1 

01 

60 

• 

83 

6 

06 

0 

00 

9 

09 

3 

03 

4 

04 

4 

04 

GTO 

22 

0 

00 

LBL 

23 

70 

D 

14 

1 

01 

• 

83 

8 

08 

5 

05 

2 

02 

GTO 

22 

0 

00 

LBL 

23 

A 

11 

80 

1 

01 

LBL 

23 

0 

00 

g 

35 

1/x 

04 

ENTER 

41 

ENTER 

41 

ENTER 

41 

X 

71 

X 

71 

90 

3 

03 

9 

09 

8 

08 

6 

06 

0 

00 

3 

03 

X 

71 

STO  1 

3301 

CL  x 

44 

6 

06 

100 

Key 

Code 

Entry 

Shown 

3 

03 

7 

07 

8 

08 

• 

83 

1 

01 

6 

06 

X 

71 

STO  4 

3304 

R/S 

84 

LBL 

23 

B 

12 

6 

06 

0 

00 

GTO 

22 

1 

01 

LBL 

23 

C 

13 

3 

03 

6 

06 

0 

00 

0 

00 

GTO 

22 

1 

01 

LBL 

23 

D 

14 

8 

08 

6 

06 

4 

04 

0 

00 

0 

00 

GTO 

22 

1 

01 

LBL 

23 

A 

11 

1 

01 

LBL 

23 

1 

01 

ENTER 

41 

X 

71 

RCL  1 

3401 

X 

71 

STO  1 

3301 

R/S 

84 

g  NOP 

3501 

g  NOP 

3501 

3501 

3501 

3501 

3501 

3501 

10 


20 


30 


40 


50 


CARD  2 


Compute  range,  h    ,  v 

r  b  max   c 


Key 

Entry 
X 
1 

CHS 

f 


Key 

Code 

Entry 

Shown 

RCL  4 

3404 

+ 

61 

STO  3 

3303 

R/S 

84 

LBL 

23 

B 

12 

STO  2 

3302 

R/S 

84 

LBL 

23 

C 

13 

STO  7 

3307 

R/S 

84 

LBL 

23 

D 

14 

g 

35 

RAD 

42 

-r 

81 

1 

01 

- 

51 

RCL  6 

3406 

■r 

81 

f-1 

32 

COS 

05 

RTN 

24 

LBL 

23 

E 

15 

g 

35 

DEG 

41 

RCL  7 

3407 

f 

31 

COS 

05 

RCL  2 

3402 

X 

71 

RCL  3 

3403 

X 

71 

ENTER 

41 

X 

71 

RCL  1 

3401 

f 

81 

STO  8 

3308 

2 

02 

RCL  3 

3403 

T 

81 

RCL  2 

3402 

ENTER 

41 

X 

71 

RCL  1 

3401 

81 

- 

51 

STO  5 

3305 

60 


70 


80 


90 


/ 
STO  6 
RCL  8 
RCL  3 

D 
RCL  8 
RCL  4 

D 

+ 

g 

ENTER 

+ 
g  x<— >y 

RCL  4 
x 

R/S 
RCL  6 

1 

+ 
RCL  5 

RCL  4 

R/S 
RCL  1 
RCL  3 


R/S 

GTO 

E 

g  NOP 

g  NOP 


IOC 


NOTE:   Before  putting  a  program  into  memory,  after  switching  to  W/PGM  mode,  press   f 
and  then  PRGM  to  clear  memory. 


CARD  3 


Compute  time-of-flight 


Key 

Code 

Entry 

Shown 

RCL  4 

3404 

+ 

61 

STO  3 

3303 

R/S 

84 

LBL 

23 

B 

12 

STO  2 

3302 

R/S 

84 

LBL 

23 

C 

13 

STO  7 

3307 

R/S 

84 

LBL 

23 

D 

14 

g 

35 

RAD 

42 

RCL  5 

3405 

X 

71 

1 

01 

- 

51 

CHS 

42 

RCL  6 

3406 

T 

81 

-1 

f 

32 

COS 

05 

RTN 

24 

LBL 

23 

E 

15 

g 

35 

DEG 

41 

RCL  7 

3407 

f 

31 

COS 

05 

RCL  2 

3402 

X 

71 

RCL  3 

3403 

X 

71 

ENTER 

41 

X 

71 

RCL  1 

3401 

81 

STO  8 

3308 

2 

02 

RCL  3 

3403 

81 

RCL  2 

3402 

ENTER 

41 

X 

71 

RCL  1 

3401 

81 

60 


70 


80 


90 


100 


Key 

Code 

Entry 

Shown 

- 

51 

STO  5 

3305 

X 

71 

1 

01 

- 

51 

CHS 

42 

f 

31 

/ 

09 

STO  6 

3306 

g 

35 

7T 

02 

ENTER 

41 

+ 

61 

RCL  4 

3404 

D 

14 

- 

51 

ENTER 

41 

f 

31 

SIN 

04 

RCL  6 

3406 

X 

71 

- 

51 

RCL  3 

3403 

D 

14 

ENTER 

41 

f 

31 

SIN 

04 

RCL  6 

3406 

X 

71 

- 

51 

- 

51 

RCL  5 

3405 

81 

RCL  1 

3401 

RCL  5 

3405 

X 

71 

f 

31 

/ 

09 

81 

R/S 

84 

GTO 

22 

E 

15 

g  NOP 

3501 

g  NOP 

3501 

g  NOP 

3501 

3501 

3501 

3501 

3501 

3501 

CARD  4 


Compute  maximizing  flight  path  angle  B 


opt 


10 


20 


30 


40 


50 


Key 

Code 

Entry 

Shown 

2 

02 

RCL  2 

3402 

ENTER 

41 

X 

71 

STO  5 

3305 

RCL  3 

3403 

X 

71 

RCL  1 

3401 

81 

- 

51 

STO  6 

3306 

RCL  4 

3404 

g 

35 

1/x 

04 

RCL  3 

3403 

g 

35 

1/x 

04 

- 

51 

RCL  5 

3405 

T 

81 

RCL  1 

3401 

X 

71 

ENTER 

41 

+ 

61 

1 

01 

+ 

61 

STO  5 

3305 

g 

35 

DEG 

41 

0 

00 

STO  7 

3307 

5 

05 

STO  8 

3308 

LBL 

23 

0 

00 

RCL  5 

3405 

RCL  3 

3403 

RCL  4 

3404 

81 

RCL  6 

3406 

X 

71 

RCL  7 

3407 

f 

31 

COS 

05 

X 

71 

gLSTx 

3500 

X 

71 

- 

51 

RCL  7 

3407 

f 

31 

60 


70 


80 


90 


100 


Key 

Code 

Entry 

Shown 

SIN 

04 

X 

71 

RCL  5 

3405 

RCL  3 

3403 

RCL  4 

3404 

81 

RCL  7 

3407 

f 

31 

COS 

05 

X 

71 

ENTER 

41 

X 

71 

- 

51 

f 

31 

/  ± 

09 
81 

1 

01 

+ 

61 

RCL  6 

3406 

a 

81 

f 

31 

/ 

09 

if"1 

32 

COS 

05 

f 

31 

TF  2 

81 

R/S 

84 

g  NOP 

3501 

RCL  7 

3407 

gx<— i>y 

3507 

STO  7 

3307 

0 

00 

g  x=y 

3523 

GTO 

22 

0 

00 

+ 

61 

+ 

61 

2 

02 

81 

STO  7 

3307 

g 

35 

DSZ 

83 

GTO 

22 

0 

00 

R/S 

84 

5 

05 

STO  8 

3308 

GTO 

22 

0 

00 

g  NOP 

',',iil 

NOTE:   Before  putting  a  program  into  memory,  after  switching  to  W/PGM  mode,  pn 
and  then  PRGM  to  clear  memory. 
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VII.   The  Appendix. 

As  indicated  in  section  HID,  the  maximum  range  for  a 
specified  cutoff  altitude   h_   and  cutoff  velocity   v_   is 
found  by  the  following  unconstrained  optimization: 


S     =  S(B   4.)  =  max  S(6) 
max     v  opt 


0<3< 


Using  the  equations  in  section  IIIA,  we  want  to  find 

6   .   such  that 
opt 


36 


=   rE  J3  (fT-fL} 


6  =  6 


=  0  . 


opt 


^   =    & 


opt 


We  find  that 


h    (VfL> 


(Y  -y  )2£-  , 

lTT  TI/ 36  ' 


where 


11 'Fm 


;r .    ~   3   r.     / 
l    2ae  \  l    / 


,  i=T,L 


For  the  derivative  of   fT~fL  with  respect  to   6   to 

vanish  at  ^    =    ^>       ,  ,  we  must  have  either: 

opt 


1. 


(YT"V6  =  6Qpt   =  0 


15 


or 


ap 

93 


=  0 


3  =  6 


opt 


We  find  that 


93 


2  2   .   00 
rTv_  sin  23 
i-i  -L 

y 


=  o 


for   3  =  0/  ±   /2,  ±  it,  etc.   Since  we  require  that  the 
launch  point  and  the  target  point  be  separated  by  the  line 
of  apses,  we  require  that  B<$<-^   .   This  requirement  is 
imposed  because  we  must  properly  interface  the  disconti- 
nuity caused  by  the  principal  angle  in  the  arc  cosine  used 
in  steps  5  and  6  of  section  III  A.   Thus   3  =  /2   is  the 
only  angle  of  interest  for  which   9P/93  =  0,  and  this 
angle  clearly  gives  rise  to  a  minimum  in  range  since  the 
trajectory  of  the  missile  is  vertical  (straight  up,  and 
then  straight  down) . 

The  value  of   3  which  maximizes   S   is  found  from  the 
requirement  that 


(VYL}3  =  3    =  °  * 
opt 


This  equation  has  defied  attempts  to  find  an  analytic 
solution,  but  it  can  be  rearranged  and  displayed  as 
equation  1  in  section  III  D,  where  the  solution  procedure 
is  detailed. 
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